T = 180e3;

t = [T*0.2,T*0.4,T*0.6,T*0.8,T];
ms = 5;

A = load('ssarun1.txt');
B = [0.003,0.038,0.050,0.063,0.110];
M = [0.149,0.093,0.059,0.031,0.012];

figure(1); clf;
plot(A(:,1), 'b'); hold on; plot(A(:,2),'r');
title('ER -> PI');
plot(t,B,'bo', 'MarkerSize', ms, 'MarkerFaceColor','b');
plot(t,M,'ro', 'MarkerSize', ms, 'MarkerFaceColor','r');

A = load('ssarun2.txt');
B=[0.184,0.101,0.066,0.023,0.015];
M=[0.005,0.075,0.078,0.088,0.273];

figure(2); clf;
plot(A(:,1), 'b'); hold on; plot(A(:,2),'r');
title('PI -> ER');
plot(t,B,'bo', 'MarkerSize', ms, 'MarkerFaceColor','b');
plot(t,M,'ro', 'MarkerSize', ms, 'MarkerFaceColor','r');